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1. Introduction 



In the lattice Landau gauge QCD simulation, we adopt two types of the gauge field definitions^] 
logf/ type: U Xtfl = et** , a\^ = -A 2# , 

U linear type: A x ^ = -{U x ^ - U x ^)\ tr i p ., 



A a 

where fr / p . implies traceless part. ( A fl (x) = i^A ll a {x)—j=, trA a A 6 = 8 ab ) 

a V2 

The corresponding optimizing functions are 
logU type: F v (g) = | |A*| | 2 = £^ * [A^A^), 
U linear type: Fu(g) = £^ tr (2 - ( + L/*/)) . 

Under infinitesimal gauge transformation g~ l 8g = £, its variation reads for either defintion as 
AFu(g) = -2(3A*|e) + <e| - dD(U g )\e) + ■■■, 
where the co variant derivativative D^(U) for two options reads commonly as 

where d^ty = 0(jc + /l) — (j>(x), and = - - -^r— _ Stationality of the optimizing function 
implies Landau gauge, the local minimum implies Gribov Region^ and the global minimum 
implies Fundamental modular(FM) region^]. 

We performed simulations on quenched configurations tabulated in Table | 



10] and un- 



quenched configurations of JLQCD[0, CP-PACS0, Columbia University [JT3[] and MILC[[T|] 
tabulated in Table 2 [fill]. 



j3 


1 /a(GeV) 


L 


aL(fm) 


definition of A 


6 


1.97 


16 


1.60 


U-linear/ log U 






24 


2.40 


U-linear/ log U 






32 


3.20 


U-linear/ log U 


6.4 


3.66 


32 


1.72 


U-linear/ log U 






48 


2.59 


U-linear/ log U 






56 


3.02 


U-linear/ log U 


6.45 


3.87 


56 


2.86 


logU 



Table 1: Configurations used in the quenched QCD simulation 



2. The Kugo-Ojima theory 



The Kugo-Ojima confinement criterion is given by the fact that the parameter c defined as 
u ab (0) = -8 ab c in the eq(0) becomes 1 



x,y 



-ip(x-y) 



(tr ( A at D 



-dD 



[A v ,A b 



(2.1) 



xy- 
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VWI 1 VWI 
am ud l am s 


N f 


l/a(GeV) 


L s 


U 


aL s (tm) 


JLQCD 


5.2 


0.1340 


0.134 


2 


2.221 


20 


48 


1.78 




5.2 


0.1355 


0.093 


2 


2.221 


20 


48 


1.78 


CP-PACS 


2.1 


0.1357 


0.087 


2 


1.834 


24 


48 


2.58 




2.1 


0.1382 


0.020 


2 


1.834 


24 


48 


2.58 


CU 


5.415 




0.025 


2 


1.140 


16 


32 


2.77 




5.7 




0.010 


2 


2.1 


16 


32 


1.50 


MILQ 


6.83(A mp ) 




0.040/0.050 


2+1 


1.64 


20 


64 


2.41 




6.76(ft mp ) 




0.007/0.050 


2+1 


1.64 


20 


64 


2.41 


MILC/ 


7.11(A mp ) 




0.0124/0.031 


2+1 


2.19 


28 


96 


2.52 








0.0062/0.031 


2+1 


2.19 


28 


96 


2.52 



Table 2: Configurations used in the unquenched QCD simulation 



The Zwanziger horizon condition[||] coincides with Kugo-Ojima criterion provided the co- 
variant derivative approaches the naive continuum limit, i.e., e/d = 1[§|. We observe that in the 
quenched simulation c saturated at about 0.8, while in the unquenched simulation it is consistent 
with 1. 





K sea or j3 


Cx 


c t 


c 


e/d 


h 


JLQCD 


K sea =0.1340 
K sea =0.1355 


0.89(9) 
1.01(22) 


0.72(4) 
0.67(5) 


0.85(11) 
0.92(24) 


0.9296(2) 
0.9340(1) 


-0.08(11) 
-0.01(24) 


CP-PACS 


K sea =0.1357 
^=0.1382 


0.86(6) 
0.89(9) 


0.76(4) 
0.72(4) 


0.84(7) 
0.85(11) 


0.9388(1) 
0.9409(1) 


-0.10(6) 
-0.05(9) 


CU 


J3 =5.415 
JS =5.7 


0.84(7) 
0.95(26) 


0.74(4) 
0.58(6) 


0.81(8) 
0.86(28) 


0.9242(3) 
0.9414(2) 


-0.11(8) 
-0.08(28) 


MILC C 


j8 =6.76 
p =6.83 


1.04(11) 
0.99(14) 


0.74(3) 
0.75(3) 


0.97(16) 
0.93(16) 


0.9325(1) 
0.9339(1) 


0.03(16) 
-0.00(16) 


MILC/ 


P =7.09 
p =7.11 


1.06(13) 
1.05(13) 


0.76(3) 
0.76(3) 


0.99(17) 
0.98(17) 


0.9409(1) 
0.9412(1) 


0.04(17) 
0.04(17) 



Table 3: The Kugo-Ojima parameter for the polarization along the spacial directions c x and that along 
the time direction c t and the average c, trace divided by the dimension e/d, horizon function deviation h 
of the unquenched Wilson fermion(JLQCD, CP-PACS), and KS fermion (MILC C ,CU,MILC/). The \ogU 
definition of the gauge field is adopted. 



3. Quark propagator 

In the unquenched lattice simulation with the improved KS fermion action, the MILC collab- 
oration has replaced the link variables by fattening^] 

Ufi (x) -► aUpix) + £ w 3 4 3 v (*) + - 
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(3) 

where is the staple contribution 



S [ ;l{x) = U v {x)U fi (x+\>)U$(x + fi) 

and added the Naik term which is a product of three link variables along one direction. 

The gauge configurations of the MILC collaboration are produced by incorporating, larger 
staple Snip an d SJlpo an d the tadpole improvement factor, and the action is called "Asqtad" action. 

In the present calculation of valence quark propagator, we do not include all improvements 
of the "Asqtad" action but incorporate staple term and the Naik term, i.e. the propagator we 
measure is the same as that of []|], and we put for Dirac operator D + mas 

1 4 (3) 

H u )x,y = -z £ T] M (x)«gn(ju)[(cif/ M (x) + w 3 £ SfivW)W 
ft=-4 v+\i 

+c 3 U^(x)Uu(x + jl)U^{x + 2jl) 8 yrX+3fl ] 
where W3 = 9/64, c\ =9/32 and C3 = —1/24, and rj^(x) is given as 



otherwise 

For KS fermion, one defines translation invariant states as 

X P ,a (x) = A=e ikx , k tl =p ll + nan 
v V 



where = Inm^/L^ (m^ = 0, (L^/2) — 1), = 0, 1. Then the Dirac operator in the tree 
level, E)(I) + m, has a proper form in the above basis as 

9 1 

{Xp>,p\®(, I )+ m \X P ,a) = ['E(%W(g sin /^ -^^Pii)+m8 aji ]8 p i p 
where the Dirac gamma matrices of KS fermions appear as 

(?n)ap = (-l) aM 4+CW,^ and 8 a p = Y[ 8 a^ 

n 

4. Calculation of the propagator 

The quark propagator is calculated by statistical average over Landau-gauge-fixed samples as 

The inversion, ^^ +m > is performed via conjugate gradient method after preconditioning as fol- 
lows. 

We define the operator ^# = (/ + —ID) with use of even- odd- sites decomposition 



m 



I ±E> 

m ° 

-V) en I 



I-L-U 
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where lower triangle and upper triangle matrices, L and U , should be properly understood. 
Using the Eisenstat trick[|]], we define 



(I-LY x Jl(I-Uy x = (I + L){I-U -L)(I + U) 
I-LU 



I 

7- (h) 2 1D P n1D nP . 



where L 2 = U 2 = understood. 

We note that ^# is hermitian and the conjugate gradient method and/or BiCGstab method are 
applicable for its inversion. Thus for solution of the equation m^tp = p , we obtain that 

= (I + U)^£-\l + L)-p 
m 

For the right hand side, p, of the equation ^#<^» = p, we put matrix site field p x = % p « (jc)/3 X 3 
where ^3x3 is the unit color source matrix. Then calculation with use of matrix site field, , {%p,a 1 ) 
yields the sample contribution to S a p (p) which is color 3x3 matrix. 

Our error estimate of the inversion calculation is 

\\^(<j> — p\\ 

n — n < a few per cent at most 

IIpII 

where the used norm is maximum norm in the space of site, color and flavor, and the accuracy gets 
10 _1 higher if L 2 norm is used. 

With use of momentum definition 

<?Ai = ^sin(/? A( )-^sin(3^ Ai ) = sin/^ ^l + ^sin 2 (^)^ , 



one obtains that 



c , x „ f ,-i{n)aSi+M{q)~8 aSi 



Since try^ = 0, trace over color and flavor yields 



tr%) = lav ffi^l ee l6N c mg) (4-1) 



On the other hand 



tr(ifqS(q)) = 16^ ^^ = \6N c q 2 ^(q) (4.2) 



The dynamical mass function of the quark is M(q) = ^ \ , and the quark wave function 

. y ^(q) 2 q 2 + SS{q) 2 
renormalization is Z2{q) = 
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5. Discussion and conclusion 

We found that the Kugo-Ojima parameter is consistent with 1 in the MILC configuration. 
The quark field renormalization Z2 is found to be infrared suppressed^, 0|. Its implication 
to the Kugo-Ojima confinement criterion is under investigation. We observed further, the renor- 



malization of Z2 and that of the running coupling are correlated QX2p- On the extraction of the 
continuum limit of the running coupling from compact lattice simulations, there is a warning from 
the Dyson-Schwinger approach[jl7|]. 
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